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Abstract 

We investigate the separability of quantum states based on covariance matrices. Separability 
criteria are presented for multipartite states. The lower bound of concurrence proposed in Phys. 
Rev. A. 75, 052320 (2007) is improved by optimizing the local orthonormal observables. 
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As one of the most striking features of quantum phenomena, quantum entanglement 
has been identified as a key non-local resource in quantum information processing such 
as quantum computation, quantum teleportation, dense coding, quantum cryptographic 
schemes, entanglement swapping and remote state preparation 1]. The study of quantum 
information processing has spurred a flurry of activities in the investigation of quantum 
entanglements. Nevertheless, despite the potential applications of quantum entangled states, 
the theory of quantum entanglement itself is far from being satisfied. One of the important 
problems in the theory of quantum entanglement is the separability: to decide whether or not 
a given quantum state is entangled. In principle the problem could be solved by calculating 
the measure of entanglement. However most proposed measures of entanglement involve 
extremizations which are difficult to handle analytically. 

There have been some (necessary) criteria for separability, the Bell inequalities 2|, PPT 
(positive partial transposition) 3| (which is also sufficient for the cases 2x2 and 2x3 bipartite 



systems |4jJ), reduction criterion 
realignment 



10 



11 



majorization criterion[7|, entanglement witnesses 



12| and generalized realignment [13 




1 as some necessary 



and sufficient operational criteria for low rank density matrices 

In [l?! by using the Bloch representation of density matrices the author has presented 
a separability criterion that is independent of PPT and realignment criteria. It is also 
generalized to multipartite case [l8]. In 19] a criterion based on local uncertainty relations 
has been presented. It has been shown that the criterion based on local uncertainty relations 
is strictly stronger than the realignment criterion 2(]|. The covariance matrices are then 



introduced to solve the separability problem in 21 



22 



23 



241 ] . Recently a criterion which 



is strictly stronger than the realignment criterion and its nonlinear entanglement witnesses 



introduced in 



201 ] has also been presented |25j. 



In this letter we study the separability problem by using the covariance matrix approach. 
An alternative separability criterion is obtained for bipartite systems. The local orthonormal 
observable dependent lower bound of concurrence proposed in [22j is optimized. The covari- 
ance matrix approach is applied to multipartite systems and a set of separability criteria is 
obtained. 

We first give a brief review of covariance matrix criterion proposed in 22j. Let Ti^ and 
71^ be d- dimensional complex vector spaces, and pab a bipartite quantum state in 7-^£g>7if . 
Let Ak (resp. Bk) be d 2 observables on TC^ (resp. 7i^) such that they form an orthonormal 



2 



normalized basis of the observable space, satisfying tr(A k Ai) = 5 k j (resp. tr(B k Bi) 
Consider the total set {M k } = {Ak <8> I <8> B k }. It can be proven that [201 ]. 



JV 2 



A' 2 



(i) 



fc=i fc=i 

The covariance matrix 7 is defined with entries 



(2) 



which has a block structure 



22]: 



7 



A C 



(3) 



where A = 7^, B = j(p B , {B k }),C i:j = (A { ®Bj) PAB - (Ai) ^(Bj) PB , p A = 

Ttb{pab), Pb = Tta(pab)- Such covariance matrix has a concavity property: for a mixed 
density matrix p = Y^PkPk with p k > and = 1, one has j(p) > J2Pkl(Pk)- 

k k k 

For a bipartite product state pab = Pa ® Pb, C in (jHJ) is zero. Generally if p^s is 
separable, then there exist states \a k )(a k \ on \bk)(bk\ on and p^ such that 



(4) 



where = ^p k ^(\a k )(a k \, {A k }), k b = ^p k j(\b k }(b k \ : {B k }). 

The so called covariance matrix criterion is made more efficient and physically plau- 
sible in 



22] . For a separable bipartite state, it has been shown that 

1 - Hp\)) + (1 - tr(p|)) 



< 



(5) 



i=l 



Criterion (jSJ) depends on the choice of the orthonormal normalized basis of observable. 

d 2 

In fact the term \Cu\ has an upper bound ||C||^ which is invariant under unitary 

i=l 

transformation and can be attained by choosing proper local orthonormal observable basis, 



where ||C||xf stands for the Ky Fan norm of C, \\C \\ kf = try CC\ with f denoting the 
transpose and conjugation. It has been shown in |26|] that if pab is separable, then 

(1 - tr(p\)) + (1 - tr(pl)) 



\\C\\ 



KF 



< 



(6) 
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From the covariance matrix approach, we can also get an alternative criterion. From ([3]) 
and (j4j) we have that if pab is separable, then 

f A- k a C \ 
X = > 0. (7) 

V C T B-K B j 

Hence all the 2x2 minor submatrices of X must be positive. Namely one has 



> 0, 



(A — ka)u Cij 

Cji (B — KB)jj 

i.e. (A — ka)h(B — Ks)jj > Cfj. Summing over all i, j and using ([I]), we get 

d? 

J2 C h - (trA-trK A ))(trB -trn B ) 

= (d - tr(p 2 A ) -d + l)(d- tr(p 2 B ) - d+ 1) = (1 - tr(p 2 A ))(l - tr(p%)). 

That is 

\\C\\h<(l-tr(p A ))(l-tr(pl)), (8) 

where ||C||hs stands for the Euclid norm of C, i.e. ||C||ijs = yJtr(CC^). 
Formulae ([6]) and jSJ) are independent and could be complement. When 

l-tr(p A Ml-tr(p B )) < \\C\\ HS < \\C\\ KF < (1 " tr(pi)) + (1 ~ Mpl)) , 

dHJ can recognize the entanglement but can not. When 

\\C\\hs < ^l-tr(p A Ml-tr(pl)) < (1 " tr(pi)) + (1 " tr(pl)) < \\C\\ KF , 

(J6]) can recognize the entanglement while (|HJ) not. 

The separability of a quantum state can also be investigated by computing the concur- 
rence. The concurrence of a pure state is given by C(|V>)) = y/2{l - Trp\) Q, where 
Pa = TrB\ip)(ip\- Let p be a state in TC^ Cg> M < N. The definition is extended to 
general mixed states p by the convex roof, 



In 



C{p) = mm <^ ^PiCiipi) : p = ^pMiipil > 

{Pit . 

\ I I ) 



(9) 



291 ] a lower bound of C (p) has been obtained, 



C(p) > \I M{M _ 1) [Max(\\T A (p)\\,\\R(p)\\) - 1] , (10) 



where Ta and R stand for partial transpose with respect to subsystem A and realignment 



respectively. This bound is further improved based on local uncertainty relations 



3l|, 



r( , . M + N-2-J2 t A%G?®I + I®G?) 

C\p) > , = , (11) 

y/2M(M-l) 

where Gf and Gf are any set of local orthonormal observables. 

Bound (fTTl) again depends on the choice of the local orthonormal observables. In the fol- 
lowing we show that this bound can be also optimized, in the sense that a local orthonormal 
observable-independent up bound of the right hand side of (fTTT) can be obtained. 

[Theorem 1] Let p be a bipartite state in Hm <S> H^. C(p) satisfies 

c( ) > 2\\C\\ KF -(l -Trp A )-(l-Trp%) 
^2M{M- 1) 

[Proof] The other orthonormal normalized basis of the local orthonormal observable space 

can be obtained from A, and Bi by unitary transformations U and V: Ai = Yl UuAi and 

i 

Bj = Y Vj m Brn- Select U and V so that C = WAV is the singular value decomposition of 

m 

C . Then the new observables can be written as Ai = UuAi, Bj = —J2 V* m B m . We have 

l m 

^Aj(l,®/ + /® J B i ) = ^[A^(l0 + A^(S i )+2((I i ®A)-(I i )( J B i })] 

i i 

= M - Trp\ + N- Trp 2 B - 2 J2( UCV ^h 

i 

= M -Trp 2 A + N -Trp 2 B -2\\C\\ KF . 



Substituting above relation to flTTTl we get (f!2l . ■ 
Bound ( fl2|) does not depend on the choice of local orthonormal observables. It can be 
easily ap plie d and realized by direct measurements in experiments. It is in accord with the 



result in [26[ where optimization of entanglement witness based on local uncertainty relation 
las been taken into account. As an example let us consider the 3x3 bound entangled state 

P=^ 9 -X>>(&l), (13) 

where J 9 is the 9 x 9 identity matrix, |f > = ^|0)(|0) - |1», |&) = ^(|0) - |1))|2), 
|6) = ^|2)(|l) - |2», |6) = ^(|1> - |2»|0), |6) = |(|0) + |1) + |2»(|0> + |1> + |2». We 
simply choose the local orthonormal observables to be the normalized generators of SU(3). 
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Bound 




FIG. 1: Lower bounds from (12) (dashed line) and (10) (solid line) 



Formula (JTOj) gives C(p) > 0.050. Formula §TB) gives C(p) > 0.052 [3l|] , while formula (JH 
yields a better lower bound C(p) > 0.0555. 



If we mix the bound entangled state (1131) with 



75 E !**)> p' = (! - x )p + 



i=0 



it is easily seen that (fl2"l) gives a better lower bound of concurrence than formula (|T0|) (Fig. 
!)• 

The separability criteria based on covariance matrix approach can be generalized to mul- 
tipartite systems. We first consider the tripartite case, Pabc £ ® ® ■ Take 
observables on Ha resp. A on Kg resp. on 7^c- Set {M^} = {Ak Cg) I ® I, / eg) A ® 
I, I ® I <g>Ck}. The covariance matrix defined by (j2J) has then the following block structure: 



( A D E^ 



7 



D T B F 
E T F T C 



(14) 



where A 



-y{p A Mk})> ^ = l(p B ,{Bk}), C 
(A^C^^-iA^^Cj) 
[Theorem 2] If pabc is fully separable, then 



\A% ) p A (-Bj ) PB , A? 



PA \ w J / PC7 J - 1 *J 



7(Pc,{C fe }), Aj = (A ® AW 

= (A ® Cj)p BC ~ (A) p B (Cj) pc- 



\D\ 


|2 

Ihs 


< 


(1 


-tr(p 2 4 ))(l 


-tr(p|)) 


\E\ 


|2 

l.ffS 


< 


(1 


"^(p 2 4 ))(l 




\F\ 


|2 

IffS 


< 


(1 


-H4))(i 





(15) 
(16) 
(17) 
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and 



2\\D\\ KF <(l-tr(p 2 A )) + (l-tr(p 2 B )), 
2\\E\\ KF <(l-tr(p 2 A )) + {l-tr(p 2 c )), 
2\\F\\ KF <(l-tr(p 2 B )) + (l-tr(p 2 c )). 



(18) 
(19) 
(20) 



[Proof] For a tripartite product state Pabc — Pa ® Pb ® Pc, D, E and F in ffT4l) are 
zero. If pabc is fully separable, then there exist states |a&}(a&| in H^, \bk)(bk\ in 7if and 
\c k )(c k \ in H%, and p fc such that 7(p) > /t4 © k.b © «c, where k a = J2pkl{\ak)(a k \, {A k }), 
k b = J2Pkl{\h)(h\,{B k }) and k c = Y,Pkl{\ck)(c k \, {C k }), i.e. 



Y 



( 



A — ka D E 
D T B-k b F 
E T F T C- k c 



\ 



> 0. 



(21) 



Thus all the 2x2 minor submatrices of Y must be positive. Selecting one with two rows 
and columns from the first two block rows and columns of Y, we have 



(A - K A )i 



Dji (B — Kb) a 



> 0, 



(22) 



i.e. (A — ka)u(B — K B )jj > \Dij\ 2 . Summing over all i, j and using (JTJ), we get 



d 2 



IDII 2 



^2 D h - ( trA ~ tr ^A)){trB - trK B ) 

(d - tr{p 2 A ) -d + l)(d- tr(p 2 B ) - d+ 1) = (1 - tr(p 2 A )){l - tr(p 2 B )), 



which proves i fTST) . ( |T6i) and (JTTj) can be similarly proved. 

From (1221) we also have (A — Ka)u + (B — K B )u > 2|Djj|. Therefore 

(trA - trK A )) + (trB - trn B ) 



< 



(d - tr(p 2 A ) -d+l) + (d- tr(p 2 B ) -d + 1) 
2 

(1 - Hp 2 A )) + (1 - tr(p 2 B )) 



(23) 



d 2 d 2 

Note that J2 D u<J2 I Ail- By using that Tr(MU) < \ \M\\ KF = TrVMW for any matrix 

i=l i=l 

0d 2 
, we have ^2 Da < \\D\\kf- 
i=i 
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Let D = WAV be the singular value decomposition of D. Make a transformation of the 

orthonormal normalized basis of the local orthonormal observable space: A, = UnAi and 

i 

Bj = ^2 V* m B m . In the new basis we have 



Dij = UuD lm V jm = (UDV% = V (24) 



lin 



Then (1231) becomes 



(l-tr(pj)) + (l-tr(p|)) 

KF S ~ 



i=i 

which proves (IT8l) . (fl9l) and (|20l) can similarly treated. ■ 
We consider now the case that pabc is bi-partite separable. 

[Theorem 3] If pabc is a bi-partite separable state with respect to the bipartite partition 
of the sub-systems A and BC (resp. AB and C; resp. AC and B), then (fl5l) . (jT6l) and (jTSI) . 
(fT9|) (resp. (HI]), ([I7D and ([IHD, (T20D; resp. (O, (JTZD and (USD, (1201)) must hold. 

[Proof] We prove the case that pabc is bi-partite separable with respect to the A system 
and BC systems partition. The other cases can be similarly treated. In this case the matrices 
D and E in the covariance matrix (fT4"l) are zero, pabc takes the form pabc = Yl.PmpA® Pbc- 

m 

Define k a = Y^Pml{PA, i A k}), k-bc = ]Cp™7(p!bc> i B k ® I J ® C k }). k bc has a form 

kbc 

where k b = J2 Pk7 {\h)(h\, {Bk}) and k c = J2Pkl(\ck)(ck\, {Cfe}), COy = J2Pm((Bi 

m 

Cj) p m c — (Bi) p m(Cj) p m). By using the concavity of covariance matrix we have 




• K A \ 



k b F 

V o (ft «o ; 

Accounting to the method used in proving Theorem 2, we get ffT5l) . ffl6|) and ffTSl . f|T9|) . ■ 
From Theorem 2 and 3 we have 

[corollary] If two of the inequalities $155, (ffSD and (ED (or ([18]), and (J20D) are violated, 
the state must be fully entangled. 

The result of Theorem 2 can be generalized to general multipartite case p G nf ®nf ® 
■ ■■®Ttf ) . Define A l a = I®/® • • • A a <g>7® • • -®J, where A = 7/ei, A Q (a = 1, 2, • • ■ d 2 - 1) 



are the normalized generators of SU(d) satisfying tr{\ a \p} = 5 a p and acting on the % 



th 



system Ti. 



d ' 



1, 2, • ■ ■ ,N. Denote {M^} the set of all A l a . Then the covariance matrix of 



p can be written as 



l(p) 



An Au ■ ■ ■ Ain 
AJ2 A22 ■ ■ ■ A2N 



\ 



A T A T 

•Sn AT /V) 



y si 1N s\ 2 n ' ' ' Ann J 

where A u = j(p, {A\}) and (^)mn = (Mn ® ^) " A><^£) for i ^ j. 

For a product state pu-N, Aij, i 7^ j, in fl25|) are zero matrices. For a fully separable 
mixed state, it has the form p = J2 Pkl^l) i^ll ® iV'fcKV'fcl ® ' ' ' ® IV^KV^I- Define 



(25) 



(26) 



Then for a fully separable multipartite state p one has 



— kau 



Aii--- 



A12 A22 - «U 2 



\ 



A T 



Ain 
A2N 



A T 



Ann — 



\ 



> 0. 



/ 



(27) 



From which we have the following separability criterion for multipartite systems: 



[Theorem 4] If a state p E H 
inequalities 



(i) 
d 



n 



(2) 



H d is fully separable, the following 



\\A 



ijWHS 



<(l-tr(^))(l-tr(^)), 



must be fulfilled for any i ^ j. 



(28) 
(29) 



We have studied the separability of quantum states by using the covariance matrix. An 
alternative separability criterion has obtained for bipartite systems, which is a supplement 
of the criterion in [22|. The covariance matrix approach has been applied to multipartite 
systems and some related separability criteria have been obtained. The local orthonormal 



observable dependent lower bound of concurrence proposed in 



22j has been optimized. 



In dealing with the multipartite cases, we have considered that all subsystems have the 
same dimensions. The results can be generalized to the case that some or all subsystems 
have different dimensions. For instance, let N max = N n be the the largest dimension. One 
can choose N% observables A^. For other subsystems with smaller dimensions, say Ni, one 
chooses Nl observables A\, k = 1, • • • , Nf, and set A\ = for k = + 1, • • • , N%. 
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